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Abstract 

A network of n wireless communication links is considered in a Rayleigh fading environment. It is assumed that 
each link can be active and transmit with a constant power P or remain silent. The objective is to maximize the 
number of active links such that each active link can transmit with a constant rate A. An upper bound is derived 
that shows the number of active links scales at most like j log n. To obtain a lower bound, a decentralized link 
activation strategy is described and analyzed. It is shown that for small values of A, the number of supported links 
by this strategy meets the upper bound; however, as A grows, this number becomes far below the upper bound. To 
shrink the gap between the upper bound and the achievability result, a modified link activation strategy is proposed 
and analyzed based on some results from random graph theory. It is shown that this modified strategy performs very 
close to the optimum. Specifically, this strategy is asymptotically almost surely optimum when A approaches oo or 
0. It turns out the optimality results are obtained in an interference-limited regime. It is demonstrated that, by proper 
selection of the algorithm parameters, the proposed scheme also allows the network to operate in a noise-limited 
regime in which the transmission rates can be adjusted by the transmission powers. The price for this flexibility is a 
decrease in the throughput scaling law by a multiplicative factor of log log n. 

Index Terms 

Wireless networks, fading channel, throughput, scaling law, random graph. 

I. Introduction 

Wireless networks consist of a number of nodes communicating over a shared wireless channel. The design and 
analysis of such configurations, even in their simplest forms, are among the most difficult problems in information 
theory. However, as the number of nodes becomes large, wireless networks become more tractable, where scaling 
laws for network parameters, such as throughput, can be derived. 

Most of the works dealing with the throughput of large wireless networks consider a channel model in which 
the signal power decays according to a distance-based attenuation law [l]-[8]. However, in a wireless environment, 
the presence of obstacles and scatterers adds some randomness to the received signal. This random behavior of 
the channel, known as fading, can drastically change the scaling laws of a network in both multihop [9]— [12] and 
single-hop scenarios [13, Chapter 8], [14]— [18]. In this paper, we follow the model of [9], [13], [19], where fading 
is assumed to be the dominant factor affecting the strength of the channels between nodes. 

In this work, we consider a single-hop network, i.e., a network in which data is transmitted directly from sources 
to their corresponding receivers without utilizing any other nodes as routers. Each communication link can be active 
and transmit with a constant power P or remain silent. Throughput and rate-per-link are the network parameters 
which are of concern to us. Despite the randomness of the channel, we are only interested in events that occur with 
high probability, i.e., with probability tending to one as n — ► oo. This deterministic approach to random wireless 
networks has been also deployed in [2], [8], [19]. 

In a previous work by the authors [19], the throughput maximization of a single-hop wireless network in a 
Rayleigh fading environment has been investigated without any rate constraints. It is shown that the maximum 
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throughput scales like log n. Also, a decentralized link activation strategy, called the threshold-based link activation 
strategy (TBLAS), is proposed that achieves this scaling law. The throughput maximization using TBLAS yields an 
average rate per active link that approaches zero as n — > oo. The same phenomenon has been observed in [1], [2], 
[9], [11]. Since most of the existing efficient channel codes are designed for moderate rates, it is a drawback for a 
system to have zero-approaching rates. Thus, from a practical point of view, it is appealing to assign constant rates 
to active communication links. In [7], it is shown that a nondecreasing rate-per-node is achievable when nodes are 
mobile. 

In this paper, we consider the problem of rate-constrained throughput maximization in a Rayleigh fading envi- 
ronment. More specifically, the objective is to maximize the number of active links such that each active link can 
transmit with a constant rate A. We derive an upper bound that shows the number of active links scales at most 
like j logn. To obtain a lower bound, first, we examine the simple TBLAS of [19] and show that it is capable of 
guaranteeing rate-per-links equal to A. The number of active links provided by this method scales like 6 (logn). 
The scaling factor is close to the optimum when A is small. However, as A grows large, the scaling factor decays 
exponentially with A, making it far below the upper bound j. This inspires developing an improved link activation 
strategy that works well for large values of desired rates, as well. To this end, we propose a double-threshold-based 
link activation strategy (DTBLAS). 

DTBLAS is attained by adding an interference management phase to TBLAS. This is done by choosing from 
good enough links only those with small enough mutual interference. The analysis of DTBLAS is more complicated 
than that of TBLAS. However, it can be carried out using some results from the random graph theory. It is shown 
that DTBLAS performs very close to the optimum. Indeed, its performance reaches the upper bound when the 
demanded rate approaches oo or 0. This shows the asymptotic optimality of DTBLAS for the rate-constrained 
throughput maximization problem. 

In all scenarios described above, the interference power is much larger than the noise power and the rates become 
independent of signal-to-noise ratio (SNR). In other words, the network performs in an interference-limited regime. 
A natural question is whether it is possible to have rate-per-links which depend on the SNR. The importance of this 
scenario, which is called the noise-limited regime, is that the transmission rate A can be adjusted by adjusting the 
transmission power P. We show that the answer to the above question is affirmative and the noise-limited regime 
can be realized by using DTBLAS. However, the throughput achieved by this method scales like l< ^°f ( ^ n , which is 
by a multiplicative factor of log logn less than what is achievable in an interference-limited regime. 

It is worth mentioning that link activation strategies studied in this paper can be considered as special power 
allocation schemes. The problem of throughput maximization via power allocation is a challenging problem for 
which only suboptimum solutions have been reported [20]-[22]. However, variations of this problem have been 
extensively studied in the literature, where the on-off scheme has frequently appeared. Recently, for a decentralized 
utility-based networlfl it is shown that the optimum power allocation follows an on-off paradigm when the number 
of links is large [23]. The on-off power allocation has been also used in [16], [17] for a cellular network in which 
the number of cell (links) are limited, but in each cell there are infinite number of users to choose from. For cellular 
systems, a distributed joint power allocation and scheduling has been proposed in [24], in which again an on-off 
strategy is followed. 

The rest of the paper is organized as follows: In Section|IIl network model and problem formulation are presented. 
An upper bound on the throughput is derived in Section [III] In Sections [IV] and |V] achievability results via 
decentralized and centralized schemes are presented. Some optimality results are provided in Section [VI] The 
operation of the network in a noise-limited regime is investigated in Section IVIII Finally, the paper is concluded 
in Section |VinI 

Notation: Af n represents the set of natural numbers less than or equal to n; log(-) is the natural logarithm 
function; [^J denotes the largest integer less than or equal to x; x 2 (M) represents the chi-squared distribution 
with M degrees of freedom; P(A) denotes the probability of event A; E(x) and Var(ir) represent the expected 
value and the variance of the random variable x, respectively; « means approximate equality; for any functions 

f(n) and h(n), h(n) — 0(f(n)) is equivalent to lim, woo \h(n)/f(n)\ < oo, h(n) = o(/(n)) is equivalent to 
lim^oo \h(n)/f(n)\ = 0, h(n) — w(/(n)) is equivalent to lim^oo \h(n)/f(n)\ = oo, h(n) = 6(/(n)) is equiv- 
alent to linin^oo \h(n) / f (n)\ = c, where < c < oo, and h(n) ~ f(n) is equivalent to rim n _, 00 h(n)/ f(n) = 1; 
an event A n holds asymptotically almost surely (a.a.s) if P(A n ) — > 1 as n — » oo. 

'Each node maximizes a locally computed network average throughput conditioned on its own channel gain. 
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II. Network Model and Problem Formulation 

The network model is the same as in [19]; however, we repeat it here for completeness. We consider a wireless 
communication network with n pairs of transmitters and receivers. These n communication links are indexed by 
the elements of Af n . Each transmitter aims to send data to its corresponding receiver in a single-hop fashion. The 
transmit power of link i is denoted by pi. It is assumed that the links follow an on-off paradigm, i.e., pi <E {0, P}, 
where P is a constant. Hence, any power allocation scheme translates to a link activation strategy (LAS). Any LAS 
yields a set of active links A, which describes the transmission powers as 

( P if IE A 

Pi = { if i£A ■ (1) 

The channel between transmitter j and receiver i is characterized by the coefficient gji. This means the received 
power from transmitter j at the receiver i equals gjiPj. We assume that the channel coefficients are independent 
identically distributed (i.i.d.) random variables drawn from an exponential pdf, i.e., f(x) = e~ x , with mean yU = 1 
and variance a 2 = 1. This channel model corresponds to a Rayleigh fading environment. We refer to the coefficients 
gu and gji (j ^ i) as direct channel coefficients and cross channel coefficients, respectively. 

We consider an additive white Gaussian noise (AWGN) with limited variance r\ at the receivers. The transmit 
SNR of the network is defined as 

P=~- (2) 
V 

The receivers are conventional linear receivers, i.e., without multiuser detection. Since the transmissions occur 
simultaneously within the same environment, the signal from each transmitter acts as interference for other links. 
Assuming Gaussian signal transmission from all links, the distribution of the interference will be Gaussian as well. 
Thus, according to the Shannon capacity formula [25], the maximum supportable rate of link i e A is obtained as 

ri(A) — log (1 + 7i(*4)) nats/channel use, (3) 

where 

7i (A) = fg (4) 

is the signal-to-interference-plus-noise ratio (SINK) of link i. 

As a measure of performance, in this paper we consider the throughput of the network, which is defined as 

T(A) = J2n(A). (5) 

i£A 

Also, the average rate per active link is defined as 

fW = (6) 

In this paper, wherever there is no ambiguity, we drop the functionality of A from the network parameters and 
simply refer to them as r^, 7,, T, or f. 

Throughout the paper, we assume all active links transmit with a same constant rate A. In this case, the 
throughput becomes proportional to the number of active links, i.e., T(A) = \A\X. Hence, the problem of throughput 
maximization becomes equivalent to maximizing the number of active links subject to a constraint on the rate of 
active links, i.e., 

max 1*41 

AQAf n 1 ' (7) 

s.t. n(A) > A, Vi e A 

This problem is referred to as the rate-constrained throughput maximization. We denote the throughput correspond- 
ing to the maximum value of this problem by T*. 

Due to the nonconvex and integral nature of the throughput maximization problem, its solution is computationally 
intensive. However, in this paper we propose and analyze LASs which lead to efficient solutions for the above 
problem. Indeed, we first show that the decentralized method of [19] is a.a.s. optimum when A is vanishingly small. 
Then, we propose a new LAS which is asymptotically optimum for large values as well as small values of A. Also, 
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for moderate values of A, there is a small gap between the performance of the proposed LAS and a derived upper 
bound. This shows the closeness of its performance to the optimum. 

For simplicity of notation, we denote the number of active links by k instead of \A\. Motivated by the result 
of [19] that shows the maximum throughput scales like \ogn, we introduce the following definitions. The scaling 
factor of the throughput is defined as 

r = lim (8) 
?woo logn 

Similarly, the scaling factor of the number of active links is defined as 

k 

k = lim . (9) 

rwoo logn 

III. Upper Bound 

In this section, we obtain an upper bound on the optimum solution of (01. This upper bound can be either 
presented as an upper bound on the throughput or as an upper bound on the number of active links. 

Theorem 1: Assume A* is the solution to the rate-constrained throughput maximization (O and k* = \A*\. Then, 
the associated throughput and the scaling factor of k* a.a.s. satisfy 

T* < log n -log log n + c, (10) 
< < ~, (ID 

for some constant c. 

Proof: For a randomly selected set of active links A with |.4| = k, the interference term Ii = ^2jeA9ji in 
the denominator of (|4]i has x 2 (2fc — 2) distribution. Hence, we have 



/•OO 

P(7* > x) = / P 
Jo 



(7i > x\h = z)f Ji (z)dz 



-xjl/p+x) z e J 

(k-2)l dZ 

e -x/p 



- ., (12) 
(1 + xj h 1 

Assume L\ is the event that there exists at least one set A C M n with |^4| = k such that the constraints in (j7]i 
are satisfied. Also, assume 70 is a quantity that satisfies A = log(l + 70). We have 



P(£i) < ( ) (P( Pi > A)) K (13) 
" * (P(7i > To))" (14) 



kj (1 + 7o)' £( ^ 1) 



(15) 



( fc ) /ne\ k e -i° k /P 

< l-r) 7^ \T7ir~TT (16) 



(1 + 70)^^^ 

_ e /c(logri-log/c-Afc+A+l-7o/p) 

where (a) is due to the union bound and (b) is the result of applying the Stirling's approximation for the factorial. 
It can be verified that there exists a constant c such that if fcA = log n — log log n + c, then, the above upper bound 
approaches zero for n — > oo. Hence, for the event C\ to have non-zero probability, we should a.a.s. have 

k A < log n — log log n + c. (18) 

This inequality holds for any feasible number of active links. By choosing k = k*, the upper bounds in the lemma 
are proved. ■ 
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IV. Lower Bound: A Decentralized Approach 

To derive a lower bound, in this section, we consider the threshold-based link activation strategy (TBLAS) 
originally proposed in [19]. 

TBLAS: For a threshold A, choose the set of active links according to the following rule 

ieA iff g u > A. (19) 

As it is seen, in TBLAS each link only needs to know its own direct channel gain. If a direct channel gain is 
above the threshold A, the corresponding link is active; otherwise, it remains silent. The value of A determines the 
achievable throughput. We show that by proper choose of the threshold, TBLAS provides a solution for the rate- 
constrained throughput maximization. The importance of TBLAS is that it can be implemented in a decentralized 
fashion. 

Let us denote the achieved throughput of TBLAS by T TBLAS . The following results are proven for TBLAS in 
[19]: 

T TBLAS ~ ne- A log(l + -^^ , (20) 

k-TBLAS ~ ne A J (21) 

\k TBLAS — ne~ A \ < £V ne~ A , a.a.s. (22) 

where the last inequality holds for any £ = u>(l). 

A necessary condition for the rate of active links being equal to A is f TBLAS = A, where f TBLAS is the average 
rate per active link achieved by TBLAS. Hence, we should choose A such that the throughput and the number of 
active links both become proportional to log n. The following lemma shows how to realize such a scenario. 

Lemma 2: Assume the activation threshold for TBLAS is chosen to be A = log n — log log n — log a for some 
a > 0. Then, a.a.s. we have 

T TBLAS = al °g( 1 + ~) (23) 

k tblas = " (24) 

f TBLAS = lOg 1 1 + - ) + 0(1). (25) 



Proof: With the specified value of A, we have ne~ A — alogn. The values of t tblas and n TBLAS are readily 
obtained by substituting this value in d20T i and ( f2TT > and using the definitions ([H} and (0, respectively. The value of 
r TBLAS is obtained by using the definition ©. ■ 

Lemma |2] indicates that by a proper choose of a, an average rate per active link equal to A is achievable; however, 
it does not guarantee that all active links can support this rate. In other words, one may ask whether TBLAS is 
capable of satisfying the constraints in problem (0. The following lemma addresses this issue and shows that a.a.s. 
the rate of all active links are highly concentrated around the average rate per active link. 

Lemma 3: Assume the activation threshold for TBLAS is chosen to be A = log n — log log n — log a for some 
a > 0. Then, a.a.s. we have 

\n-f\ < 2.^^(1 + 0(1)), Vi e A, (26) 
y or log n 

where f = log ( 1 + — V 
V a J 

To prove the lemma, we need the following result about the central limit theorem (CLT) for large deviations. 

Theorem 4 ( [26]): Let {Y m } be a sequence of i.i.d. random variables. Suppose that Y\ has zero mean and finite 
positive variance v and satisfies Cramer's condition^. For Z m = ^= Y^jLi define F m (y) — P(Z m < y). If 
V > 0, y — 0(m 1/6 ), then 

1 - F m (y) = [1 - *(„)] cxp f-^£=) + O [ j , (27) 

2 A random variable Y satisfies the Cramer's condition if its moment-generating function exists in some interval with the center at the origin. 
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where $(y) is the cdf of normal distribution and 63 = E(Y 1 3 ). 

Proof of Lemma [3} From the definition of r.i and the concavity of the log(-) function, we have 



< 



Thus, to prove the lemma, it is sufficient to prove that a.a.s. 



1 




7» 


a 





or equivalently 
where 



' log log n 
a 3 log n 

X- < Ji < x + 



(l + o(l)), Vi e A, 



x± =- 1 ±21 



' log log n 
alogn 



(l + o(l)) • 



Here, we just prove the left-side inequality in (|3T1 i. The other side can be proved in a similar manner. 
Let £2 denote the event that 

7i > X-, Vi e A. 

In the following, we show that ¥(£2) — > 1 as n — > 00. 

Denoting the cdf of 7$ conditioned on |„4| = k by F 7 (x, k), the probability of the event £2 is obtained as 



P(£ 2 



(a) 



(6) 



^P(|^|=fc)P(£ 2 ||^|=fc) 

n 



fc=0 



> P(\A\ = k)(l-F J (x^,k)) k 



k=k- 



(c) 



> (l-F 7 0r_,fc + )f + ^ P(W = A) 



(l-F^fc+^+P^ < |^| <fc+), 



(28) 
(29) 

(30) 
(31) 
(32) 

(33) 

(34) 
(35) 
(36) 

(37) 
(38) 



where (a) is because the channel gains are independent, (b) is valid for any < fc_ < k + < n and (c) is due to 
the fact that (1 — F y (x, k)) k is a decreasing function of k. According to (l22l i. by choosing 



-A 



a 



log n ± £ y a logn, 



(39) 
(40) 



for some £ — > 00, we have P(fc_ < \A\ < k + ) — > 1. Hence, to prove P(£2) ~~ * 1> it is enough to show that 
(1 — F 1 (x-, k + )) k+ — > 1. However, due to the inequality 

N&4 



it is enough to show that 



(1 - F 7 (x_, k+)) + > 1 - fc+F 7 (x_, fc+), 
fc_l__F 7 (a;_, fc + ) — > 0. 



(41) 
(42) 
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To prove (|42| >. we provide an upper bound on k + F 1 (x~ 1 k + ) and show that it approaches zero as n — ► oo. We have 



F„(a;_, k + ) 



(a) 



P(7i<a-IW 
/ 

P 



9i; 



V 

/ k 4 



(b) 
< 



E 



9j- 



> 



> 



(43) 



/ 



where (a) is based on A = {1, ■ • ■ , fc + }, which has been assumed for simplicity of notation, and (b) is due to the 
fact that, in TBLAS, gu > A for any i e A. Let us define Yj — — 1, which has the variance v = 1. Thus, 



the right-hand-side (RHS) of (l43l translates to the complementary cdf of Z = 



1 



rewritten as 



where 



1 



Eyii^. i-e. 63) can be 



K r (x-,k+) < l-P(Z<y), 



£-i-(*4--i) 



By substituting A = log n — log log n — log a and the value of x_ from d32| ) into j45l ), we obtain 



y = 2Vloglogn(l + o(l)). 



(44) 



(45) 



(46) 



Since Yj is a shifted exponential random variable, its moment-generating function exists around zero and the 
Cramer's condition is satisfied. Also, by choosing m = fc+ — 1 we have y — (^(m 1 / 6 ). Hence, the result of Theorem 
|4]can be applied to calculate the complementary cdf of Z. Consequently, by using d27l i with #3 = E(Y?) = 2, d44l) 
can be rewritten as ... 



F y (x-,k+) < [l-$(y)]exp - 



Noting that y 3 = o(^/k+) and using the inequality 1 — $(y) < — - 



O 



(47) 



from and (1461 , we conclude that 



y 



(48) 

= exp(-loglogn(l + o(l))). (49) 

It is clear that the above upper bound approaches zero as n — > 00. Hence, P(£2) 1 and the proof is complete. 

■ 

Lemma|3]shows that with the specified threshold for TBLAS, all active links can transmit with rate A = log(l+— ). 
Hence, TBLAS provides a solution, albeit suboptimum, for the problem (|7). Lemmas [2] and [3] reveal the following 
relation between the demanded rate A and n TBLAS as well as t tblas 

(50) 
(51) 



e A - 1 



e A - 1 



Noting that for small values of A, the RHS of d50b can be approximated as j and using the upper bound in 
Theorem [U it turns out that TBLAS is close to the optimum for small values of A. 
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V. Lower Bound: A Centralized Approach 

Although TBLAS enjoys the simplicity of decentralized implementation, its performance is far from the optimum. 
This can be seen by comparing the upper bound in Theorem[T]and the achievability result in ( T50b - A reason for this 
suboptimality is that the mutual interference of the active links is not considered in choosing A. In this section, we 
provide an LAS that performs close to the upper bound in Theorem [T] and turns out to be asymptotically optimum 
when A is very large or very small. We name this method double-threshold-based LAS (DTBLAS). 

DTBLAS: For the thresholds A and 5 

i. Choose the largest set A\ C J\f n such that ga > A for all i € Ai. 

ii. Choose the largest set A2 Q A\ such that < 8 and < 6 for all i, j E A 2 . 
The set of active links is A = A 2 . 

This strategy chooses the links to be active in a two-phase selection process; in the first phase, which is basically 
similar to TBLAS, a subset Ai of the links with good enough direct channel coefficients is chosen. In the second 
phase, which is the interference management phase, a subset of links in Ai is chosen such that their mutual 
interferences are small enough. Note that the second phase of the strategy requires full knowledge of the channel 
coefficients. Hence, this scheme should be implemented in a centralized fashion. 

We aim to find A and S such that the throughput is maximized subject to the rate constraints of the active links. 

For simplicity, we use the notation fej = \Ai\ for i = 1, 2. Without loss of generality, assume A4 = {1, • • • , fc^}. 
By using (O, ©, and (O, and applying the Jensen's inequality, the throughput is lower bounded as 

T > k 2 log (l + 1 A fca ) , (52) 

where Z; = Y^!j=i 9ji- Since < 5, the mean and variance of Ii depend on S. More precisely, we have 

E(Ji) = (fe - (53) 
Var(Ji) = (fc 2 -l)«r 2 , (54) 

Se~ s 

A = I'- {.'/;, .'/,;. < 5} = 1 - - _ (55) 

a 2 = Vax{g ji \g j i<5} = l- S _ 6 e _ s (56) 

Assume 6 is a constant and k% — > 00 as n — > 00. To simplify the RHS of d52b . we apply the Chebyshev inequality 
to obtain the upper bound 

-. k 2 

-rJ2 Ti < (*a-i)A + V>, (57) 



where 



-2 . . 
1=1 

which holds a.a.s. for any ip — w(l). Consequently, the lower bound (l52l becomes 

T > k 2 log (1 + — — ) a.a.s. (58) 

Note that the constant 1/p — p, is absorbed in the function ip. Since ip can be chosen arbitrarily small, say with an 
order smaller than pk 2 , we can rewrite ( 1581 as 

T > T (59) 

W ^ 6re / A 

^dtb las = fc 2 (log (1 + — J + o(l) J a.a.s. (60) 



denotes the throughput achievable by DTBLAS. 

Since is a random variable, the right hand side of d60l l is a random variable as well. However, the following 
discussion shows that k 2 is highly concentrated around a certain value. Hence, it can be treated as a deterministic 
value. 
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Construct an undirected graph G(Ai, E) with vertex set A\ and the adjacency matrix E = [e^] defined as 



1 ! 9ij < S and gji < 6 
; otherwise 



The probability of having an edge between vertices i and j, when gji and g^ have exponential distribution, equals 



P={l-e~ 5 y. (61) 

The definition of G implies that G E G(ki,p), where Q(ki,p), which is a well-studied object in the literature [27], 
is the family of k\ -vertex random graphs with edge probability p. 

In the second phase of DTBLAS, we are interested to choose the maximum number of links whose cross channel 
coefficients are smaller than 8. This is equivalent to choosing the largest complete subg rapr0 of G. The size of the 
largest complete subgraph of G is called its clique number and denoted by cl(G). The above discussion yields 

fc 2 =cl(G), for some G&Q(ki,p). (62) 

Although the clique number of a random graph G is a random variable, the following result from random graph 
theory states that it is concentrated in a certain interval. 

Theorem 5: Let < p < 1 and e > be fixed. The clique number cl(G) of G S G(m,p), for large values of 
m, a.a.s. satisfies si < cl(G) < s 2 where 

Si = L21og b m-21og b log b 7n(l-p) + 21og b (e/2) + l + (-l) l e /pJ, i = 1, 2, (63) 

b = l/p. 

Proof: The theorem is a direct result of Theorem 7.1 in [28], which states a similar result for the stability 
number of random graphs. Using the fact that the stability number of a random graph G(m, p) is the same as the 

clique number of a random graph G(m, 1 — p), the theorem is proved. ■ 

Corollary 6: Consider DTBLAS with parameters A and S. The number of active links, k DTBLAS = a.a.s. 
satisfies k'_ < k DTBLAS < k' + , where 

k' ± = L21og b ne- A -21og b log b n e - A (l-i) + 21og b (e/2) + l±e/p + (l)J (64) 
and b = (1 - er s )- 2 . 

Proof: According to d22l . a.a.s. we have ki = ne~ A + 0(^V ne~ A ). Assuming £ = o(V ne~ A ), and by 
substituting this value of k\ into d62l and using Theorem |3J the corollary is proved. ■ 

The next lemma indicates how to choose the thresholds A and S such that the throughput and the number of 
active links both become proportional to log n. As a result, a constant average rate per active link is achieved. 

Lemma 7: Assume the threshold A for DTBLAS is chosen to be 

A = (l-a')logn(l + o(l)), (65) 
for some a' > and 6 is a constant. Then, a.a.s. we have 

-a' 



2 



"DTBLAS 



DTBLAS 



log(l-e- 5 )' 

-rv' 



log(l 



( (l-a')log(l- e - A -; ^ 



(66) 



(67) 



= log(l-(^«^] +0 (l). (68) 



a'p, 

3 A complete graph is a graph in which every pair of vertices are connected by an edge 
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Proof: For the number of active links, we have 

(«) 

(6) -a'(l + o(l)) 



k DTBLAS ~ 21og b 7ie A (69) 

logn, (70) 



log(l-e- 5 ) 

where (a) is based on Corollary [6] and (b) is obtained by using d65l ). From ( f70b . and by using the definition (|9), 
k dtblas ls obtained as given in d66l l. 

The number of active links in ( T70b can be used along with the value of A in (l65l l to rewrite ( f60b as 



T 

DTBLAS 



log(l — e -5 ) \ a'/i 



logn. (71) 



The scaling factor t dtblas , as given in the Lemma, is obtained by using the value of /t from d55l ) and applying 
the definition (|8). The value of r DTBLAS is obtained by using the definition (|6j. This completes the proof. ■ 

According to this lemma, by proper choose of the constants a' and S, the average rate per active link r DTBLAS 
can be adjusted to be equal to the required rate A. A natural question is whether, under the specified conditions in 
DTBLAS, all active links can support the rate A. The following lemma addresses this issue and shows that a.a.s. 
the rate of all active links are highly concentrated around the average value f DTBLAS . 

Lemma 8: Consider DTBLAS with thresholds S and A = (1 — a') logn for some a' > 0. Then, a.a.s. we have 



for some constant c > 0, where 



/logTogn 
V log™ 

= log A_(L-«')log(l-e-^ 



a fj, 

Proof: See Appendix lAl ■ 
According to Lemmas [7] and [8] when maximizing the throughput of DTBLAS, S should be a constant and A is 
obtained from another constant a'. Hence, the rate-constrained throughput maximization simplifies to an optimization 
problem with constant parameters a' and 6. Assume 70 is a quantity that satisfies A = log(l + 70), i.e., 70 is the 
required SINR by the active links. Instead of the number of active links, we can maximize the scaling factor of 
the number of active links given in Lemma [7] Hence, the rate-constrained throughput maximization problem (|7]) is 
converted for DTBLAS to the following optimization problem 

max - - —rr (73) 



,8 log(l-e~ A ) 

(l-a')log_(l_-e- 5 ) 
Se- 



tt- / e —s X = 7o- (74) 



1 

1 - e" 

Note that in contrast to problem (|7), there is only one constraint in this problem. However, according to Lemma [8] 
this single constraint guarantees the required rate for all active links. From the equality constraint ( f74b . parameter 
a' can be found in terms of 5 as 

of -log(l-e-*) 

Tb^l-^TpyJ-loga-e-') 

By substituting this value in the objective function d73l . we obtain the following equivalent unconstrained 
optimization problem 

min 70 - Y^—s^j - log (l - e' 5 ) . (76) 

Consequently, (a 1 * , 5*), the solution of (f73l . can be obtained by first finding S* from d76l > and then substituting it 
into d75l l to obtain a'* . Due to the complicated form of d76i >. it is not possible to find S* analytically and it should 
be found numerically. 
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Fig. 1. Optimum of the threshold 8 and the parameter a' vs. the demanded rate A. 



Fig. Q] shows S* and a'* versus A. The values of 5* and a'* can be replaced in d67l i and d66l l to obtain the 
maximum throughput scaling factor (j* DTBLAS ) as well as the maximum scaling factor for the number of active 
links ( K * DTBLAS )- The value t* dtblas is shown in Fig. [2] Depicted in the figure is also the throughput scaling factor 
of TBLAS obtained from (l5TT l. As it is observed, for small values of A, the performance of TBLAS and DTBLAS 
are almost the same. However, as A grows larger, the scaling factor of TBLAS approaches zero, but the scaling 
factor of DTBLAS approaches 1. This shows some kind of optimality for DTBLAS which will be later proven 
formally. Figure [3] demonstrates the tradeoff between the number of supported links and the demanded rate-per-link 
for TBLAS and DTBLAS. The tradeoff curve for TBLAS is obtained from (ISO) . The upper bound from Theorem 
Q] is also plotted for comparison. As observed, for a ceratin value of A, DTBLAS can support larger number of 
users, especially for larger values of A. Indeed, the tradeoff curve of DTBLAS is very close to the upper bound. 
Specifically, for large values of A, these two curves coincide. This will be later proven formally. 



VI. Optimality Results 

Although the behaviour of DTBLAS is numerically described in Figs. [HE] and[3] it is possible and also insightful 
to obtain closed form expressions for 5* and a'* as well as n* DTBLAS and t* dtblas when A is very small or very 
large. An important result of these extreme-case analyses is the asymptotic optimality of DTBLAS. 

Setting the derivative of the objective function d76l > equal to zero reveals that, at the optimum point, 6, satisfies 

e A (l-e _a -6) + 6 = 0. (77) 

Two extreme cases of large A and small A are discussed separately in the following. 
a) Large A: In this case, solving ( ITTT i yields 

5* = 2e~ A + O (e" 2A ) . (78) 
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Consequently, a'*, t*. and k* are obtained as 



i-**m + o('\ (so) 



A \ A 2 

-x + °{^)- (81) 

As it is seen from the above equations, for large values of A, 6* becomes very small and a'* approaches one. This 
means, when large rate-per-links are demanded, it is more crucial to manage the interference than to choose links 
with high direct gain. 

b) Small X: In this case, solving ( l77b yields 

5* = j + \ + 0(A). (82) 
Consequently, a'*, t* dtblas , and k* dtblas are obtained as 

a'* = e-W (l + l + (\)\ (83) 
r* DTBLAS = l-^ + 0(A 2 ) (84) 

< TBLAS = + ^ 

The above equations show that for small values of A, S* is very large and a'* is very small. In other words, 
DTBLAS is converted to its special case, TBLAS. 

The above discussion yields the following optimality result on DTBLAS. 

Theorem 9: Consider the rate-constrained throughput maximization problem ||7). Assume r* and k* are the 
maximum achievable scaling factors of the throughput and the number of supported links, respectively. Also, 
assume r* and k* are the maximum scaling factor of the throughput and the number of active links 

DTBLAS DTBLAS b b ^ 

when DTBLAS is deployed. Then, a.a.s. we have 

(86) 
(87) 
(88) 

A ' ' ' 

lilll DT B LAS = J (89) 

Proof: The proof of the theorem is straightforward by using the upper bounds provided in Theorem Q] and the 
asymptotic achievability results provided in this section. ■ 



(r* 

V DTBLAS 


r*) 


-o, 


\ K DTBLAS 


o 


= o, 


(t* 

\ DTBLAS 


<) 


= o, 


K* 

DTBLAS 




— 1 



VII. Noise-Limited Regime 

In the previous sections, we considered an interference-limited regime in which the noise power is negligible in 
comparison with the interference power. In this case, the achievable throughput is not a function of the network 
SNR. In other words, changing the transmission powers does not affect the supportable rate of each link. However, 
in a practical scenario, it is appealing to have rates which scale by increasing p. This way, the transmission rates 
can be easily adjusted by changing the transmission powers. Specifically, it is desirable that the rate of active links 
a.a.s. scale as 

ri = Iog(l + i7 ^-), Vie A, (90) 

for some /3j = 0(1), which are the design parameters. At the same time, we require the conditions of problem (Q, 
i.e. Ti > A, be satisfied. In this section, we show how to realize such a situation by using DTBLAS. 
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According to d90l >, we should a.a.s. have Ii = j3i, where Ii is the interference observed by active link i and is 
defined in ( T52l . However, this requires that E(ij) = /3j. Noting that E(ij) = (k-2 — l)/t (see d53"l>). we conclude that 
all /3jS should take a same value, say (3. Hence, a necessary condition for being in the noise-limited regime is 

(fca - l)/t = /3, (91) 

where /? = 0(1) is a design parameters. Later, we show that d9TT > is also a sufficient condition for operating in a 
noise-limited regime. 

Based on the above discussion, we propose the following scheme for choosing the parameters of DTBLAS for 
a noise-limited regime: For a given required rate A = log(l + 70) and the interference /3, 

i. choose A as 

A = A = 7o (l/p + /3). (92) 

ii. choose 6 such that d9"Tl ) is satisfied. 

Note that the selection of A is such that the rate constraints T"j > A are satisfied. Also, as will be shown later, the 
selection of 8 is such that operation in the noise-limited regime is guaranteed. 

The next step is to solve d9"Tl > to obtain the value of S and the corresponding number of active links fca- By using 
d55l >. which gives the value of a in terms of S, it is clear that d9"TT l holds only if 6 — ► as ki — * 00. In this case, 

d55l > converts to p, = — + 0(8 2 ) and d9lT > simplifies to 

k 2 S = 2/3 a.a.s (93) 

To solve d93l and obtain S, we should first obtain the value of k 2 in terms of n and 6. From d22i > and condition 
(1921 . the number of links chosen by phase (i) of DTBLAS is obtained as 

fci = ne~ A ° + O UVruF^} . (94) 

Also, recall from d62l that k^ is the clique number of a random graph G(ki, p), where p is obtained from d6lT ). 
Since 8 — > 0, d6TT i can be rewritten as 

p = S 2 + 0{6 3 ), (95) 

which approaches zero as well. Note that Theorem|5] which was adopted from [28], and a similar result that appears 
in [29], are valid only for a fixed value of p. A natural question is whether a similar concentration result on the 
clique number of random graphs holds when p approaches zero. In the following lemma, we address this issue and 
obtain a concentration result on the clique number for zero-approaching values of p. 

Lemma 10: Let p = p(m) be such that p = o(l) and p = Lu{m~ a ) for all a > 0. For fixed e > 0, the clique 
number cl(G) of G £ G(m,p) a.a.s. satisfies [s\ < cl(G) < |_sj + 1, where 

s = 2 log,, m - 2 log 6 log b m + 1 - 41og b 2 - - — -, 

logo 

b = l/p. 

Proof: See the Appendix. ■ 
By using this lemma, d94l l, d95l ), and assuming £ = o(y / n), the number of active links a.a.s. becomes 



k 2 = 

Thus, ( |93l can be rewritten as 



log n — log log n 



-log 5 
log n — log log n 



(96) 



-log S 

Assuming | log /3| = o(loglogn), it can be verified that the solution of d97l > is 

2/3 log log n 



■5 = 2(3. (97) 



log n 

With this value of S, the number of active links is obtained from d96l ) as 

log log n 



(l + o(l)). (98) 



(99) 
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As mentioned before, we should show that the selected values of S and A for DTBLAS, yields the network to 
operate in the noise-limited regime. The following theorem addresses this issue. 

Theorem 11: For the values of A and 5 given in (|92l i and d98t , respectively, the interference of active links a.a.s. 
satisfy 

\Ii - /3\ -> 0, Vi e A. (100) 
Proof: By using the central limit theorem it can be shown that 

\h-p\< — . , Vie A, (101) 

Vlogn 

which readily yields the desired result. Since the calculations are similar to those in the proof of Lemmas [3] and [8] 
we omit them for brevity. ■ 
Lemma 12: Let T NL denote the throughput achieved by DTBLAS in the noise-limited regime described above. 
Then, almost surely we have 

Jos f 1 + tAt?) < lim < lo s f 1 + • ( 102 > 

Proof: According to Theorem QT| the throughput is obtained as 

^=|>g(l + T7 ^). (103) 



>fc 2 log( 1+ - A ° ) ■ (104) 



Due to the fact that ga > A , we have 

T > k, lncr I 1 J , 

The left-hand-side inequality in the lemma is readily obtained by using this inequality and the value of &2 from 
d99i >. For the right-hand-side inequality, by utilizing the Jensen's inequality in d 1 03b , we obtain 

T - sw °4 + ^§?iP)- <105) 

According to the law of large numbers and due to the fact that gu > A , we have 

T Yl 9u -> E(g l4 |3. i2 > A ) = 1 + A . (106) 

%— 1 



The result is obtained by using dl05l ), ( 11061 ), and the value of k<z from 

It is observed that the price for operating in the noise-limited regime is a decrease in the throughput by a 
multiplicative factor of log log n. 

VIII. Conclusion 

In this paper, wireless networks in Rayleigh fading environments are studied in terms of their achievable 
throughput. It is assumed that each link is either active and transmits with power P and rate A, or remains silent. 
The objective is to maximize the network throughput or equivalently the number of active links. First, an upper 
bound is derived that shows the throughput and the number of active links scale at most like log n and j log n, 
respectively. To obtain lower bounds, we propose two LASs (TBLAS and DTBLAS) and prove that both of them 
a.a.s. yield feasible solutions for the throughput maximization problem. In TBLAS, the activeness of each link is 
solely determined by the quality of its direct channel. TBLAS, which can be implemented in a decentralized fashion, 
performs very close to the upper bound for small values of A. However, its performance falls below the upper bound 
when A grows large. In DTBLAS, the mutual interference of the links are also taken into account when choosing 
the active links. It is demonstrated that DTBLAS not only performs close to the upper bound for A — > 0, but its 
performance meets the upper bound when A — > oo. The above discussions take place in an interference-limited 
regime in which the transmission power P does not affect the transmission rate A. However, we show that by a 
proper choose of the DTBLAS parameters, the rate-constrained network can also operate in a noise-limited regime; 
this feature of the DTBLAS comes at the price of decreasing the network throughput by a multiplicative factor of 

log log 71. 
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Appendix A 
Proof of Lemma[8] 

The proof is based on the same arguments as in the proof of Lemma|3] Thus, here we just highlight the differences. 
Let us define 7 as 

(l-qQlogfl 
a fj, 

Similar to the proof of Lemma [3] it is enough to show that a.a.s. 

x'_<-yi<x' + , (108) 

where 

with d = c/7. We only prove the left side inequality in d 1 081 >; the other inequality can be proved in a similar 
manner. 

Let £3 denote the event that 

%>x'_, Vie A, (110) 

In the following, we show that P(£s) — > 1 for some d > 0. 

Note that with A = (1 — a') logn, the parameter k' + in Corollary [6] is obtained as 

k' + = n DTBLAS logn- aloglogn (111) 

< k otblas logn, (112) 

where K OTBLAS is given in d66l l and a > is a constant. Denoting the cdf of 7$ conditioned on |„4| = k by 
Fj(x, k), we have 

P(As) > (l-F^x'_,k' + )f+J>(k'_ < \A\ <k' + ) (113) 

« (l~F 7 (x'_,k' + )f+ (114) 

> (l~F 7 (x'_,n DTBLAS \ogn)y—— losn , (115) 

where (a) is obtained in the same manner as d38l . (b) results from Corollary [6] and (c) is due to ( II 121 > and the fact 
that (1 — F~f(x, k)) k is a decreasing function of k. To show that the RHS of d 1 15b tends to one, we upper bound 
k dtblas lognF 7 (j;L,K D3 , BtAS logn) and show that it approaches zero. 
Similar to the derivation of j43l , it can be shown that 



F 1 (x'_, k dtblas logn) < P 



-DTBLAS l0 S™ 



x A 1 



(116) 



Let us define Yj = gji — fi, where fi is obtained from (|53T >. Random variable Yj has the variance v = a 2 , where 
a 2 is given in d56l l. By defining Z = — ^^^^^=^^^^^= Y^-^ BLAS log " Yj, (II 16b can be reformulated as 

V^DTBLAS^Sn- 1) J^i 

Ey(xr_,K BTBLA8 1ogn) <l-P(Z<y), (117) 

where 

p~ - \ ~ Kdtblas logn - 

y = — — . — . (118) 

\/( K DTBLAS log"" 1)<7 2 

By substituting A = (1 — a') logn and the value of x'_ from ( 11091 ) into ( II 18b , we obtain 

= jJ k^sP (l + O ( * )V (119) 

V (t z \ \ y' log n log log n / / 
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It is straightforward to show that the moment-generating function of Yj exists around zero. Hence, the Cramer's 
condition is satisfied. Also, by choosing m = n DTBLAS logn — 1, the condition y = C^to 1 / 6 ) is satisfied, as well. 
As a result, Theorem [4] can be utilized to calculate the RHS dll7b as 

1 - ?{Z < y) = [1 - exp ( d3y3 =) + O ( —L=L= ) (120) 

V W V K DTBLAS J \ V K OTBLAS l0 S " / 



By combining ( II 17b . (1120b . and ( 11 19b , and noting that 63 is a constant, y 3 — o(y / K DTBLAS logn), and 1 — <P(y) < 

-y 2 /2 

, we conclude that 



e 2 

k dtblas lognFj(x, k dtblas logn) < k dtblas logn—— (121) 

= exp ( (1 - ijlv^MlL ) log log n + 0(log log logn) 
\ zcr z 

It is clear that if d is chosen large enough, the above upper bound approaches zero as n — > 00. This completes the 
proof. 

Appendix B 
Proof of Lemma[To1 

The proof is based on the standard second moment method. 

A. Preliminary Calculations 

Assume Y s is the number of cliques of size s in G. Let us denote its mean and variance by u s and a 2 s , respectively. 
According to [29], we have 



s /s\ fm—s\ 

^W^Z (6 G)_i), (123) 

£=2 V s / 

where 6 = 1/p. By applying the Stirling's approximation to d 1 22b , we obtain 



n s i s\ { m—s\ 

U? 



m 1 2 3(3-1) 
"s = T rP^~ (124) 

V&+*(m-s) m - s+ 3 



- rzr (I- ir p ^ (125) 

For any e > 0, the inequality 1 — x > e~( 1+e ^ x holds for sufficiently small values of x. Since we are interested in 
small values of s/m, from this inequality and ( 1124b . we obtain 

U s < e s( - l °S m-log s+(l+e)- log fc) (126) 

Equation ( 1123b is readily converted to the following inequality 

where ^ _ 

F <? = — — & w . (128) 
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By using the definition of the binomial coefficients, we obtain 



Fe < 2*- (m ~ s)! . (m ~ 5)! S —. b ^ (129) 

< t J -r s b— (130) 

(m — s) 8 

/TO \ — l(t-l) 

= 2 s ■ ( 1) (131) 



s 

Noting that — 3> 1, the above inequality can be approximately written as 

^<2 s -(-) -6—. (132) 

\ TO/ 

Using ( I127l i and ( 1132b . we obtain 

4<E^), (133) 

where 

S (*) = s log 2 + £(log s - log to + - log - - log b) (134) 

is a quadratic convex function with a minimum at £q — ~ Toff §• define 

s o = 2 log b to - 2 log b log b to - 2 log b 2. (135) 
It is easy to show that if s > sq, then g(s) > g(2). Hence, (113 31 > can be simplified as 

2 

^| < e log S +g( S )_ (136) 
Ms 

According to the Markov's inequality, we have 

P{F s >l}<Ma- (137) 

For a fixed e > 0, define 

si = 21og b TO-21og b log b TO + l + 21og 6 (e/2) + -^-. (138) 

logo 

Using ( 1126b . it is easy to verify that for s > si, we have fi s — > as m — > oo. Hence, from (1137l i. we conclude that 

P{Y S >1}^0, for s>si (139) 

as to — ► oo. This means a.a.s. the clique number of G is less than si, i.e., we have the following upper bound on 
cl(G) 

cl(G) < si a.a.s. (140) 
According to the Chebyshev's inequality, we have 

2 

P{n = 0}<^§. (141) 

K 

For a fixed e > 0, define 

s 2 = 21og b TO-21og b log b TO + l-41og b 2- (142) 

log 

Using ( 1136t , it is easy to verify that for s < S2, we have of//*, — * as to — > oo. Hence, from (114U . we conclude 
that 

P{F s =0}^0, for s<s 2 (143) 
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as 771 — * oo. This means a.a.s. the clique number of G is not less than [s 2 J, i.e., we have the following lower bound 
on cl(G) 

cl(G) > [s 2 \ a.a.s. (144) 

For sufficiently small e, the difference between the upper bound s% and the lower bound s 2 is less than one. Hence, 
from (|1401 > and (11441 ) we can conclude that 

Ls 2 J < cl(G) < Ls 2 J + 1 a.a.s. (145) 
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